Fluorine-induced local magnetic moment in graphene: A hybrid DFT study 
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Recent experimental evidence that fluorinated graphene creates local magnetic moments around F adatoms 
has not been supported by semilocal density-functional theory (DFT) calculations where the adsorption of an 
isolated F adatom induces no magnetic moment in graphene. Here, we show that such an incorrect prediction 
of the nonmagnetic ground state is due to the self-interaction error inherent in semilocal exchange-correlation 
functional. The present hybrid DFT calculation for an isolated F adatom on graphene predicts not only a 
spin-polarized ground state with a spin moment of ~1 ^b^ but also a long-range spin polarization caused by 
the bipartite nature of the graphene lattice as well as the induced spin polarization of the graphene states. The 
results provide support for the experimental observations of local magnetic moments in fluorinated graphene. 

PACS numbers: 71.15.Mb, 73.20.Hb, 73.22.Pr, 75.75.-c 



Control of the electronic and magnetic properties of 
graphene is of significant technological interest in emerging 
spintronics applications such as new types of magnetic sen- 
sors as well as spin-based random access memory. ^'^ In partic- 
uler, the adsorption of adatoms on graphene has been utilized 
to change the unique electronic band structure of graphene, 
offering an effective way to create local magnetic moments 
in nonmagnetic graphene. ^~^ Such a chemical functionaliza- 
tion of graphene with adatoms can thus be useful to manipu- 
late spin transport in graphene, which is essential to the op- 
eration of spintronics devices. For this respect, there have 
been a number of theoretical and experimental studies to ex- 
plore adatom-induced magnetism in graphene ^^ However, 
the precise identification and characterization of such mag- 
netism in graphene are still at an early stage and sometimes 
remain controversial. For instance, experimental and theoreti- 
cal studies of fluorinated graphene have reported inconsistent 
results J^~— 

Recently, Hong et al^ observed a colossal negative 
magnetoresistance in dilute fluorinated graphene, suggesting 
adatom-induced magnetism. Subsequently, they observed an 
anomalous saturation of the phase coherence length in the 
same system, which provided evidence for spin-flip scatter- 
ing caused by an adatom-induced local magnetic moment. ^^ 
By magnetization measurements as a function of externally 
applied magnetic field, Nair et al^ observed spin- 1/2 para- 
magnetism in the fluorinated graphene (CF^) samples with x 
increasing from 0.1 to 1.0. Since the measured number of 
paramagnetic centers is three orders of magnitude less than 
the measured number of F adatoms on graphene, ^^ they in- 
terpreted the paramagnetic centers in terms of local magnetic 
moments arising from the edges of F clusters which are eas- 
ily formed on grapheneJ^^ii Therefore, these recent experi- 
mentsi^r— evidenced the presence of adatom-induced local 
magnetic moments in fluorinated graphene with a wide range 
of F concentration. 

In contrast to the above-mentioned experimentaii^rii evi- 
dences of F-induced local magnetic moment in graphene, pre- 
vious DFT studies with semilocal treatment of exchange and 
correlation, i.e., generalized gradient approximation (GGA), 
predicted that an isolated F adatom on graphene has a non- 



magnetic ground state.--— According to scanning tunneling 
microscopy experiments, a single-atom defect on graphene 
gives rise to a quasilocalized state at the Fermi level, which 
extends over several nanometers around the defect.i24i^ We 
note that semilocal GGA functionals tend to stabilize ar- 
tificially delocalized electronic states due to their inherent 
self-interaction error (SIE),^^~^i since delocalization reduces 
the spurious self-repulsion of electron itself. In this sense, 
previousi^"— GGA calculations for an isolated F adatom 
on graphene seem to overestimate the relative stability of 
a rather delocalized non-magnetic configuration compared 
to a spin-polarized configuration, resulting in a nonmag- 
netic ground state. Therefore, it is interesting to examine 
if the observedi^r— local magnetic moment in fluorinated 
graphene can be predicted by the correction of the SIE with 
an exchange-correlation functional beyond the GGA. 

In this Rapid Communication, we present a new theoret- 
ical study for an isolated F adatom on graphene, which ex- 
tends the previous work by considering a hybrid exchange- 
correlation functional. Unlike previous GGA calculations pre- 
dicting a nonmagnetic ground stater^"— the present hybrid 
DFT calculation predicts that an isolated F adatom induces 
a spin moment of ~1 ^b in graphene. Especially, we find 
that F adsorption gives rise to a long-range spin polarization 
of Pz orbitals in graphene as well as a band-gap opening of 
graphene 71 bands. Our prediction of a spin-polarized ground 
state is consistent with recent experimental ^^~^^ observations 
of F-induced local magnetic moments in graphene. 

The present first-principles DFT calculations were per- 
formed using the FHI-aims code^^ for an accurate, all-electron 
description based on numeric atom-centered orbitals, with 
"tight" computational settings and accurate tier-2 basis sets. 
For the exchange-correlation energy, we employed both the 
GGA functional of Perdew-Burke-Ernzerhof (PBE)^^ and the 
hybrid functional of Heyd-Scuseria-Emzerhof (HSE).^^ For 
simulation of the adatom-graphene system, we used a periodic 
sheet geometry in which one adatom adsorbs on one side of 
graphene within a 3\/3x3\/37?30° unit cell (see Fig. 1) and 
adjacent graphene sheets are separated with a vacuum spac- 
ing of 50 A. In order to remove the dipole interactions be- 
tween periodic sheets, we applied corrections to the electro- 
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FIG. 1 : (Color online) Top view of the optimized structure of the F- 
graphene system. Large (bright) and small (dark) circles represent C 
and F atoms, respectively. The inset highlights the F-adsorbed region 
with the side view. The unit cell and the Voronoi cell are shown as 
the solid and dashed lines, respectively. The numbers in C atoms 
represent the atomic site for the magnetic moment m/ plotted in Fig. 
3(c). 



static potential and the total energyi2^ For the Brillouin zone 
integration, we used 12x 12 k points in the surface Brillouin 
zone. All the atoms were allowed to relax along the calcu- 
lated forces until all the residual force components are less 
than 0.02 eV/A. 

We begin to optimize the atomic structure of an isolated F 
adatom on graphene using the PEE functional. Here, we use 
our calculated graphene lattice constant of 2.467 A, which is 
very close to the experimental^^ value of 2.46 A. The opti- 
mized structure of such an F-graphene system is displayed in 
Fig. 1. It is seen that the C atom (labeled as Co) bonded 
to the F adatom shows a protrusion from the graphene layer, 
due to the change of bonding character from sp2 to sp^,- 
like hybridization. The vertical displacement of Co from the 
graphene layer and the bond length Jf-Cq (<^Ci-C2) between 
the F (Co) and Co (Ci) atoms are calculated to be 0.510 and 
1.553 (1.482) A, in good agreement with previousi^"— PBE 
calculations. Figure 2(a) shows the calculated band structure 
for the F-graphene system. We find that the Fermi level £'f is 
shifted lower in energy relative to the Dirac point of graphene 
by 0.45 eV and the F-induced state lying at the Fermi level is 
partially occupied, indicating some charge transfer from the 
graphene % states to the F-induced state. Indeed, our analysis 
of the Mulliken charges shows a charge transfer of 0.30 e from 
graphene to the F adatom, representing a partial ionic bonding 
between the F adatom and graphene. We note that the charge 
character [see the inset of Fig. 2(a)] of the F-induced state 
represents a quasilocalized electron which is distributed dom- 
inantly over the sites of the sublattice (labeled (3) opposite to 
the one (labeled a) to which an isolated F adatom is attached. 

To support the recent experimentaU^ril evidences that flu- 
orinated graphene exhibits local magnetic moments, we per- 
form the additional spin-polarized calculation for an isolated 




FIG. 2: (Color online) Spin-unpolarized band structures of the (a) F- 
graphene and (b) H-graphene systems computed using the PBE func- 
tional. The spin-polarized band structure of the F-graphene system 
computed using the HSE functional is given in (c). The band disper- 
sions are plotted along the symmetry lines of the Brillouin zone of 
the 3\/3x3\/3/?30° unit cell [see the inset in (a)]. The inset in (a) 
also shows the charge character of the F-induced state at the K point 
with an isosurface of 0.002 el^ , together with a perspective view 
around the F adatom. The energy zero represents the Fermi level. 
The majority and minority bands in (c) are drawn with the solid and 
dashed lines, respectively. 



F adatom on graphene using the PBE functional. However, 
we are unable to find the stabilization of a spin-polarized con- 
figuration, consistent with previous GGA calculations J^~— 
On the other hand, it is noticeable that our PBE calcula- 
tion for an isolated H adatom on graphene (employing the 
same 3>/3x3>/37?30° unit cell as in the calculation of the F- 
graphene system) predicts a spin-polarized ground state with a 
magnetic moment of 1 ^b per adatom, in good agreement with 
previous GGA calculations. -"i^il^ Figure 2(b) shows the band 
structure of the spin-unpolarized solution for the H-graphene 
system. It is seen that the H-induced state has a very narrow 
band width of 0.06 eV, leading to a splitting of the spin-up and 
spin-down bands (see Fig. 1 of the Supplemental Material^i) 
via the Stoner-type ferromagnetic exchange^'^^ as discussed 
below. In contrast, the F-induced state has a relatively larger 
band width of 0.61 eV [see Fig. 2(a)], thereby giving rise to 
the spin-unpolarized ground state. 

The fact that the F-induced state is largely delocalized over 



C atoms of the p sublattice leads us to speculate that the SIE 
inherent to the PBE functional may cause the incorrect pre- 
diction of the spin-unpolarized solution as a ground state. We 
note that the local or semi-local DFT scheme tends to overes- 
timate the stability of delocalized states due to the SIE:— ~i^ 
that is, the LDA or GGA functional within which an elec- 
tron interacts with its own charge density unphysically lowers 
the energy of delocalized states, thereby incorrectly predict 
some materials to be metals rather than insulators. In order to 
correct the SIE, we use the hybrid HSE functional^^ to per- 
form the spin-unpolarized and spin-polarized calculations for 
an isolated F adatom on graphene. Our hybrid-DFT calcu- 
lations predict a spin-polarized ground state with a magnetic 
moment of 0.97 jub per adatom, which energetically favors 
over the spin-unpolarized configuration by 32 meV. There- 
fore, we can say that the HSE functional corrects the spurious 
delocalization to reproduce the observed^^"^^ local magnetic 
moment in fluorinated graphene. 

The spin-polarized band structure of the F-graphene sys- 
tem obtained using the HSE calculation is displayed in Fig. 
2(c). We find that (i) the F-induced state is almost fully spin- 
polarized to yield a magnetic moment of ~1 /lib per adatom, 
(ii) the TT states of graphene have a band gap opening of 
0.51 eV at the F point, and (iii) the graphene states are spin- 
polarized with a spin splitting of less than 0.45 eV. Here, the 
spin polarization of the graphene states is attributed to their 
exchange interactions with the spin moment of the F-induced 
state. Consequently, the F-graphene system has a long-range 
spin polarization [Fig. 3(a)] not only due to the bipartite na- 
ture of the graphene lattice but also to the exchange spin- 
polarization effect in the valence states. Figure 3(c) shows 
that the calculated magnetic moment mt within the muffin-tin 
sphere centered at each C atom in the a ((3) sublattice gives 
a negative (positive) spin polarization. The values of mt are 
-0.01, -0.03, -0.02, and -0.02 (0.07, 0.03, and 0.02) ^5 for 
Co, C2, C4, and C^ (Ci, C3, and C5), respectively. Note that 
the magnetic moment within the muffin-tin sphere^^ at the F 
adatom shows a positive spin polarization with m = 0.06 jub. 

In Fig. 3(c), we also plot the values of m^ in the H-graphene 
system obtained using the HSE and PBE calculations. We find 
that the HSE values of m/ in the H-graphene system are similar 
to those in the F-graphene system, indicating that the magni- 
tudes of exchange splitting in the H- and F-induced states are 
close to each other. It is notable that, for the H-graphene sys- 
tem, the PBE values of m^ are relatively smaller than the corre- 
sponding HSE ones, caused by a relatively smaller exchange 
splitting of the H-induced state in the former compared to in 
the latter as discussed below. As a consequence, PBE gives 
rise to a relatively shorter range of spin polarization [see Fig. 
3(b)] compared to the HSE result (Fig. 2 in the Supplemental 
Material). 

Figure 4(a) shows the density of states (DOS) for the spin- 
unpolarized configuration of the F-graphene system obtained 
using the PBE functional, together with the result of the H- 
graphene system. We find that the F-graphene system has the 
DOS ofD = 2.14 states/eV per spin at £^f, which is smaller 
than the corresponding one (D = 4.00 states/me V) for the H- 
graphene system. This reduction of D in the F-graphene sys- 
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FIG. 3: (Color online) Spin densities of the (a) F-graphene and (b) 
H-graphene systems obtained using the HSE and PBE functional, 
respectively. The majority (minority) spin density is displayed in 
dark (bright) color with an isosurface of 0.02(— 0.02) e/A^. The inset 
in (a) enlarges the F-adsorbed region with the side view. For the F- 
graphene and H-graphene systems, the calculated magnetic moment 
rui within the muffin- tin sphere centered at the C/ atom (see Fig. 1) 
is plotted in (c). Here, the muffin-tin sphere radius is taken as 0.71 



tem does not allow for the presence of a ferromagnetic insta- 
bility.^^ On the other hand, the H-graphene system satisfies 
the Stoner criterion for the ferromagnetic instability: i.e., DI 
> 1, where the Stoner parameter / can be estimated as ~0.26 
eV by the spin splitting of the H-induced state (see Fig. 4(b) 
and Fig. 1 in the Supplemental Material^) obtained using 
PBE. Such a spin splitting of the H-induced state increases 
to ~0.86 eV by the HSE calculation, as shown in Fig. 4(b). 
This increase of exchange splitting in HSE gives the relatively 
larger values of mt compared to the PBE values [see Fig. 3(c)] . 
The resulting HSE value of / = 0.86 eV in the H-graphene sys- 
tem is slightly larger than the corresponding one (/ = 0.77 eV) 
in the F-graphene system, which is estimated by the spin split- 
ting of the F-induced state [see Fig. 2(c) and Fig. 4(b)]. Thus, 
we can say that the HSE calculation for the F-graphene sys- 
tem fulfills the Stoner criterion with DI = 1.99 > 1 to yield a 
ferromagnetic ground state. ^^ 

The present HSE result that an isolated F adatom induces 
a magnetic moment of 0.97 jub in graphene can be considered 
in reasonable agreement with the experimental observation of 
magnetic F adatoms with spin 1/2.^^ It is, however, noticeable 
that the experimentii did not detect a sign of ferromagnetic 
ordering but observed a purely paramagnetic behavior. Here, 
only one out of ~ 1 ,000 F adatoms was estimated to contribute 
to the paramagnetism. This observed paramagnetism was ex- 
plained by clustering of F adatoms, such that neighbouring 



(a) 



> 

S 

en 
O 

Q 



F-graphene(PBE) 
\ H-graphene(PBE) 




2.5 



-2.5 

-5 



F-graphene(HSE) 
/ \ / \ H-graphene(HSE) 
( \ H-graphene(PBE) 




-J \ \ \ \ \ \ \ L 

1 

Energy (eV) 



FIG. 4: (Color online) Calculated density of states for the (a) spin- 
unpolarized and (b) spin-polarized configurations of the F-graphene 
and H-graphene systems obtained using the PBE or HSE functional. 
The energy is relative to the Fermi level. The Gaussian broadening 
scheme with a width of 0. 1 eV is used. 



adatoms residing on different graphene sublattices do not con- 
tribute to the overall magnetic moment. We note that the av- 
erage spacing between paramagnetic centers located at cluster 



edges was estimated to be ^10 nm,^"^ which is larger than the 
distance (1.3 nm) between neighboring F adatoms in our em- 
ployed supercell. Thus, in order to measure a ferromagnetic 
ordering in the F-graphene samples, it would be required to 
reduce the spacing between magnetic moments either by de- 
creasing the cluster sizes which may be controlled by the sizes 
of ripples in graphene ^^ or by adsorbing F adatoms at a suf- 
ficiently low temperature to avoid the tendency towards clus- 
tering. 

In summary, using a first-principles DFT calculation with 
the HSE exchange-correlation functional, we investigated the 
spin polarization for an isolated F adatom on graphene, which 
has not been adequately described by previous semilocal DFT 
calculations. ^^~^^ The present hybrid DFT calculation showed 
that an isolated F adatom induces a spin moment of ~1 jub 
in graphene, contrasting with the semilocal DFT prediction 
of a nonmagnetic ground state. We found that F adsorption 
gives rise to a long-range spin polarization of p^ orbitals in 
graphene due to the bipartite nature of the graphene lattice and 
the exchange spin-polarization effect in the graphene states. 
Thus, our prediction of the spin-polarized ground state in the 
F-graphene system provides theoretical support for the recent 
experimental observations^^""— of local magnetic moments in 
fluorinated graphene. 
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